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A novel quantum ground state of matter is realized in a bosonic model on three dimensional fee 
lattice with emergent low energy excitations. The novel phase obtained is a stable gapless boson 
liquid phase, with algebraic boson density correlations. The stability of this phase is protected 
against the instanton effect and superfluidity by self-duality and large gauge symmetries on both 
sides of the duality. The gapless collective excitations of this phase closely resemble the graviton, 
although they have a soft lo ~ k 2 dispersion relation. There are three branches of gapless excitations 
in this phase, one of which is gapless scalar trace mode, the other two have the same polarization 
and gauge symmetries as the gravitons. The dynamics of this novel phase is described by a new set 
of Maxwell's equations. The defects carrying gauge charges can drive the system into the superfluid 
order when the defects are condensed; also the topological defects are coupled to the dual gauge field 
in the same manner as the charge defects couple to the original gauge field, after the condensation of 
the topological defects, the system is driven into the Mott Insulator phase. In the 2 dimensional case, 
the gapless soft graviton as well as the algebraic liquid phase are destroyed by the vertex operators in 
the dual theory, and the stripe order is most likely to take place close to the 2 dimensional quantum 
critical point at which the vertex operators are tuned to zero. 
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I. INDRODUCTION 

Ground states of quantum many-body systems can be- 
have qualitatively differently from classical states. Classi- 
cal ground states like ferromagnetic order and Neel order 
can survive from quantum fluctuation at zero temper- 
ature. The quantum ferromagnetic state and quantum 
Neel ordered state can be described in a similar way as 
their classical counterparts. Searching for nonclassical 
spin states can be traced back to the early proposal of 
the RVB state on frustrated lattices^, which was shown 
to be of great importance to the High Tc superconduc- 
tivity in cuperates. 

All the classical orders break certain symmetry, either 
internal symmetry or space symmetry. Quantum spin 
ground states without classical orders are termed spin 
liquids. In order to make sure the quantum ground state 
has no tendency to order, the low energy emergent gauge 
symmetry is usually applied, since as is well known, the 
gauge symmetry cannot be broken without condensation 
of matter fields^. After two decades of work, several types 
of nonclassical ground states have been identified. For 
instance, the spin liquid with Zi gauge symmetry has 
been realized in either quantum dimer model on triangu- 
lar lattice^, or spin- 1/2 model on Kagome lattice 5 . The 
existence of the Zi spin liquids in 2+1 dimensional space 
is based on the fact that the quantum Zi gauge theory 
has a deconfined phase, which manifests itself as a disor- 
dered spin state with fractionalized spin excitations. In 
these models, the ground states have no classical order, 
i.e. there is no symmetry breaking. The ground state 
degeneracy depends on the topology of the space mani- 
fold, which is termed topological order. The excitations 
contain the deconfined gapped Zi gauge charges and the 



Zi vortices (which are usually called visons). 

Another type of spin liquids contain gapless collective 
excitations (which are different from the gapless Magnon 
excitations in spin ordered state), and the spin-spin cor- 
relation functions fall off algebraically, i.e. the state is in 
a stable critical phase. These algebraic spin liquids usu- 
ally involve U(l) gauge field, or U(l) gauge field interact- 
ing with gapless matter fields. The stability of algebraic 
spin liquids is very tricky, because algebraic spin liquids 
are critical, presumably there are supposed to be relevant 
perturbations which can drive the system into an ordered 
phase. Actually in most cases the algebraic liquid phases 
are fine-tuned, for instance, the Rokhsar-Kivelson (RK) 
point of quantum dimer model on square lattice^. Due 
to the presence and proliferation of monopoles in 2+1 di- 
mensional systems^, the gapless photon excitations are 
generally gapped out, and the matter fields are confined. 
In the 3+1 dimensional space, the compact gauge the- 
ory has a deconfined photon phase. Based on this result, 
several photon liquid phase has been proposedi2*ii*i£. In 
the 2+1 dimensional systems, the monopoles are almost 
always proliferating. The algebraic liquid phase is very 
hard to survive. The monopoles are only irrelevant at 
certain critical point, for instance the RK point of the 
quantum dimer models, as well as the transition between 
valence bond solid phase and Neel state for spin- 1/2 an- 
tiferromagnetic systemi&ii. Recently it has been argued 
that at the large N level (N is the number of flavors of 
matter fields), there is a stable algebraic liquid phase in 
2+1 dimensional spaced. 

Gapless bosonic excitation is one property which char- 
acterizes criticality. Since here the critical phase is stable, 
it implies that even without any continuous symmetry 
breaking there is a gapless bosonic excitation invulner- 
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able to perturbations. This is actually a quite amazing 
property. As is well-known to all, the Goldstone the- 
orem is one way to protect the gaplessness of bosonic 
systems^. By breaking continuous global symmetry, the 
coset of the unbroken symmetry subgroup corresponds 
to the gapless modes, which are called Goldstone modes. 
Almost all the stable gapless bosonic excitations are re- 
lated to certain continuous symmetry breaking, for in- 
stance, the phonon in solids is related to the spatial 
translational symmetry breaking, the magnon is related 
to the SU(2) spin symmetry breaking. However, the gap- 
less critical excitations of algebraic spin liquid phase do 
not rely on any symmetry breaking. Therefore the alge- 
braic spin liquid phase kindled the questions like " original 
lights". It has been proposed that the photon, which is 
one of the most fundamental particles in the universe, 
might not be so fundamental, it could be collective exci- 
tations of lattice spin modelsAi. 

All the stable algebraic liquid phases which have been 
found so far involve U(l) gauge field theory. This U(l) 
gauge symmetry does not exist in the high energy (micro- 
scopic) model, instead it emerges at low energy Hamilto- 
nian, due to the effective constraint V • E = p imposed 
by the spin interaction, p is the background static charge 
distribution. Although the microscopic models which 
have been proposed so fa»i2iii*i£ are different on the lat- 
tice scale, the 3+1 dimensional photon liquid phases are 
all the same at long scale. In one previous piece of work 17 
and the current work, a new type of algebraic spin (bo- 
son) liquid phase has been realized. Based on the stan- 
dard spin-boson mapping (S z = n—h, and S* ~ exp(z#)), 
the model is presented in the bosonic version. The alge- 
braic boson liquid state studied in this paper is a new 
state of matter, it broadens the family of algebraic liq- 
uid phases. The new algebraic liquid phase does not 
involve U(l) gauge theory. Instead, the gauge symme- 
try of this model is identical to the gauge symmetry of 
linearized Einstein gravity. There are three branches of 
gapless collective excitations in this new bosonic alge- 
braic liquid phase. One of the collective excitations is 
a scalar mode, the other two gapless collective excita- 
tions have the same gauge symmetry and polarizations 
as the gravitons. However, the graviton excitations in 
our theory have a softened quadratic dispersion, uj ~ k 2 . 
The spin-spin (boson density) correlation functions fall 
off algebraically, with an exponent bigger than that of 
the 3+1 dimensional photon liquid. Like gravitational 
theory, the basic variables in this theory are symmet- 
ric rank-2 tensor, which in the gravity language are the 
linearized metric tensor. The dynamics of the graviton 
phase is described by a new set of Maxwell's equations, 
with the vectors E and B replaced by rank two symmet- 
ric tensors. The charges and topological defects enter 
the Maxwell's equations in a special form, and the con- 
densations of charges and topological defects will drive 
the liquid phase into the superfluid phase and the Mott 
insulator phase respectively. 

This paper is organized as follows. In the second sec- 



FIG. 1: The flipping term of the dimer model. 

tion, a brief review of the 3+1 dimensional photon spin 
liquid phase is presented. The discussion of the gravi- 
ton spin liquid will follow the same logic as the photon 
spin liquid phase. In the third section, the 2+1 dimen- 
sional version of the graviton model is discussed. In the 
2+1 dimensional case, the graviton theory is dual to a 
scalar boson model with quadratic dispersion. However, 
the vertex operators will generally gap out the graviton 
excitations. In the forth section, the 3+1 dimensional 
graviton model is described, and it is shown that the 
graviton phase (Gaussian phase) is self-dual, and hence 
stable. In section V, more properties of this model are 
derived. The boson density correlation function is calcu- 
lated, and a new set of Maxwell's equations are supposed 
to describe this Gaussian phase. In section VI, we dis- 
cuss the possible experimental realization of this model, 
as well as the future work. If the model is written in terms 
of fermionic operators instead of bosonis operators, novel 
non-fermi liquid behavior is expected. It is also possible 
to develop this theory to be a new candidate of quantum 
gravity theory. 



II. REVIEW OF 3+1 DIMENSIONAL PHOTON 
LIQUID PHASE 

Let us first briefly review some basic properties of the 
3+1 dimensional photon liquid phase as a warmup. Al- 
though many different models have been proposed in the 
last few yearsi2*ii*ia, the phases obtained have very sim- 
ilar properties at long scale. Therefore, let us take the 3 
dimensional Quantum Dimer Model on cubic lattice as an 
example. On the cubic lattice, every site is shared by six 
links, and only one of those links is occupied, by exactly 
one dimer. On every square face, if two parallel links 
are occupied by dimers, they can be resonated to dimers 
perpendicular to the original ones (Fig. QJ. Besides this 
resonating term, another diagonal weight term for each 
flippablc plaquette is also included in the Hamiltonian 

# = E H>(= I + h - c -) + V (\MW + !=><= D- « 

The dimer model can be mapped onto a rotor model. A 
rotor number can be defined on each link to describe the 
presence or absence of dimers: n = 1 when the link is oc- 
cupied by a dimer, and n = when the link is empty. The 
Hilbert space of this quantum system is a constrained 
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FIG. 2: The sign convention of mapping from rotor number 
onto the vector electric field E in XY plane. 

one, with the constraint $^i=i n% = \ around each site. 
Let us define quantity E^a as E^ = {—l) lx+iy nifi (The 
sign distribution is shown in Fig. [5J . Now, the constraint 
of this system can be rewritten as the Gauss's law for elec- 
tric fields, diEi = ±1. Notice that because the quantity 
Ei_a is defined on links, a natural vector notation can be 
applied: E^ = (E^, -E^g, E^ Here the derivatives are 
all defined on the lattice diEj(x) = Ej(x + i) — Ej(x). 

The background charge distribution ±1 on the cubic 
lattice plays a very important role in the solid phase, i.e. 
the confined phase. The form of the crystalline phase 
can be determined from the Berry phase induced by the 
background charge distributioni2*Ai. However, as we are 
focusing on the algebraic liquid phase, the background 
charge is not very important. Let us instead impose the 
constraint diEi = 0. Because of this local constraint, 
the low energy physics will be invariant under the gauge 
transformation A — > A + V/ (A is the conjugate variable 
of E, which is related to the phase angle 8 associated with 
each rotor number n through Ai(x,i) = (—^-Y m+ly 9 X {), 
which is exactly the gauge symmetry of the U(l) gauge 
theory. Now the effective Hamiltonian of this theory 
should be invariant under the gauge transformation. The 
Hamiltonian can be effectively written as 

H = — icos(V x A) + l/(2n)E 2 (2) 

This Hamiltonian is the 3 dimensional compact QED, 
which should have a deconfined photon phase 8,9 . In this 
phase, the cosine functions in J2J can be expanded, and 
the system can be described in the following Gaussian 
fixed point Lagrangian 

L = (d T A) 2 — c 2 (V x A) 2 . (3) 

The constraint \J -E = introduces a Lagrange multiplier 
A to the Lagrangian A (Vi ■ Ei). After integrating out 
E, the gauge invariance can be enlarged to 3+1 dimen- 
sional space time, A^ — > A^ + d^f, with /i = 0,1,2,3. 



The 3+1 dimensional Lagrangian reads 

L ~ -F^F^. (4) 

The effective Lagrangian for the photon phase has been 
derived above, and the reason for the existence of the 
photon phase is actually the remarkable self-duality of 
this photon phase and the gauge symmetry, as discussed 
below. 

First, because the constraint diEi = is strictly im- 
posed on this system, the matter field, i.e. the de- 
fect which violates this constraint, is absent. As is well 
known, the f (1) gauge symmetry can be spontaneously 
broken due to the condensation of the matter fields. How- 
ever, without matter field, the local gauge symmetry can- 
not be broken spontaneously^. The superfluid phase of 
the original system implies (9) ^ 0, written in the low en- 
ergy variables, it reads (A) ^ 0. The nonzero expectation 
of vector potential A breaks the local gauge symmetry. 
Therefore the dimer superfluid order is ruled out. 

The other possible instability of the photon phase is 
towards the gapped solid phase, which can be analyzed 
in the dual theory. We can introduce the dual vector h 
and dual momentum vector 7? (h and 7? are both defined 
on the faces of the cubic lattice) as 

E = V x h, V x A = 7T. (5) 

. One can check the commutation relation and see that 
h and 7? are a pair of conjugate variables. The Gauss's 
law constraint on this system is automatically solved by 
the vector h. Now, the field theory for the photon phase 
is self-dual: 

L = {d T A) 2 - c 2 (V x A) 2 , V • E = 0, 
L = (d T h) 2 -c 2 (V xh) 2 ,V -tt = Q. (6) 

The violation of the dimer constraint (for instance the 
hole of the doped dimer model) can be viewed as the 
gauge charges. The defects couple to the gauge field vec- 
tor in a gauge invariant manner 

L e = -tcos{S79 (e) -!) + ■■■. (7) 

#( e ) is the phase angle of the defect creation operator, 
which plays the role of the electric charge in the QED 
language. There are usually two flavors of matter fields, 
since besides the U(l) gauge symmetry, there is an extra 
global U(l) symmetry, which corresponds to the global 
conservation of total holon number. When the defects 
condense, the gauge bosons are gapped out by Higgs 
mechanism, the system enters the superfluid order, and 
the global U(l) symmetry is spontaneously broken and 
becomes the gapless Goldstone phason mode in the su- 
perfluid phas o 18 ' 19 . 

In principle, a vertex operator cos(2itNh) is supposed 
to exist in the dual Hamiltonian, due to the fact that E 
only takes on integer values. N is an integer depending 
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on the Berry phase of the vertex operator. When this 
vertex operator is relevant, it will gap out the photon 
excitation and drive the system into a crystalline phase, 
according to the Berry phase. However, the vertex op- 
erator is irrelevant in the photon phase. Notice that, 
because the theory is self-dual, the dual theory has the 
same gauge invariance h —> h + V/ as the original the- 
ory. Also, vector 7? is subject to the same constraint as E, 
diiri = jnj. However, the vertex operator cos(2iv Nh) 
breaks the gauge symmetry and thus the correlation func- 
tion between two vertex operators is zero at the Gaussian 
fixed point, i.e. the vertex operator is irrelevant in this 
Gaussian theory. 

Because the theory is self-dual at the Gaussian phase, 
the magnetic monopoles (the dual charges) should couple 
to the dual vector potential h in the same manner as 
the coupling between electric charge and original vector 
potential A 

L m = -t' cos(2vr NV9 {m) ~ 2wNh). (8) 

N is again introduced by the Berry phase, corresponding 
to the multi-monopole event. Unless the dual charges 
(the monopole) condense and break the dual U(l) gauge 
symmetry, the photon phase is always stable. The gap- 
lessness of this phase is protected by both the self-duality 
and the gauge symmetry. 

The photon phase is an algebraic liquid phase, the 
dimer density-density correlation function falls off alge- 
braically. When t = V in equation QJ, the lattice model 
can be solved exactly, and the equal-time density corre- 
lation falls off asiS 

(n(O)n(r)) ~ 1/r 3 . (9) 



III. GRAVITON MODEL IN 2 DIMENSIONAL 
SPACE 

Let us begin with the 2 dimensional example. This 
2 dimensional model is built on a square lattice with 
quantities defined on both sites and centers of plaquettes. 
Let us assume there is one orbital level on the center of 
each plaquette, and two orbital levels on each site (see 
FigEJ. The Hamiltonian of the system contains 3 terms, 
H = Hq + Hi + H%. Hi is merely the nearest neigh- 
bor hopping term between the sites and the centers of 
plaquettes 

2 

#i = E £ — *( & I.* 6 s + h - c -)- ( 10 ) 

<i,j> a=1 

Here i denotes the site of the lattice, and j denotes the 
center of plaquette, the summation is over all the hop- 
pings between each site and its 4 nearest neighbor pla- 
quettes. H 2 is an on-site interaction H 2 = U(n — n) 2 , 
which fixes the average filling per orbital state, n is the 
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FIG. 3: The structure of the 2d lattice, on each site there are 
two orbital levels, the occupation number is (711,712). On each 
face center there is one orbital level, with occupation number 
n. 

average particle number per site, for simplicity it is taken 
to be 1. 

The most important term in this Hamiltonian is Hq. 
It is a two body interaction between particle numbers 
with a special form. Each link of the square lattice is 
shared by two plaquettes and two sites. We denote the 
link between sites i and i + x as (i,x), and denote the two 
plaquettes connecting to the link &s i = i + l/2x + l/2y 
and i — y — i — l/2y + l/2x. The term in Hq which 
involves this link reads 

V(n i+ i/2x+i/2$ + n t -i/ 2 y+i/2x + 2ni,j + 2n 1A+x ~ 6h) 2 , 

(11) 

and the interaction term in Hq involving the link (i, y) is 

V(n i+1 /2x+y + n l _ 1 / 2x+ i/2y + 2n 2 .i + 2n 2 ,i +y - 6rl) 2 . 

(12) 

Notice that, for links in x direction only ri\ is in this 
interaction term, and for links in y direction only n 2 is 
involved. If V is much bigger than t, the summation in 
the brackets in both Ijlll) and l|12l) should be zero, this 
becomes a constraint on the Hilbert space if V is large. 

Because the square lattice is a bipartite lattice, we can 
stagger the sign for each sublattice. Let us define new 
variables E xx (i) = r)i(rii t i - 1), E yy (i) = r/ 4 (n 2 ,i - 1); and 
E xy (i) = ru{ n i ~ !)■ V s are signs defined on sites and 
centers of plaquettes, rji(rji) = ±1, with + for sublat- 
tice A (A) and — for sublattice B (B). The constraint 
approximately imposed by Hq can now be rewritten as 

2Q X E XX -f- QyE X y 0, 

9 X E X y + 2dyEyy = 0. (13) 

Again all the derivatives are defined on lattice. The con- 
vention of the staggered signs is depicted in Fig. 

The nearest neighbor hopping t term will generate cer- 
tain "ring exchange" term by perturbation theory, which 
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FIG. 4: The sign convention in the definition of symmetric 
tensor Eij. After introducing the signs on the lattice, the 
constraint imposed by Ho can be written compactly as 1131 . 

is allowed by the constraint (|llf> and i|12|) at low en- 
ergy. However, without doing the perturbation literally 
one can guess the form of the ring exchange term from 
the form of the constraint (|13|l . Just like the constraint 
diEi = generates gauge transformation Ai — > Ai + di<p, 
the current constraint on (|13() will generate gauge 
transformation for its conjugate variable A^ 

Aij^Aij + difj + djfi, (14) 

and the low energy ring exchange term generated from 
perturbation theory should be invariant under this gauge 
transformation. Aij is related to the phase angles of the 
original boson creation and annihilation operators by in- 
troducing staggered sign distribution r\i. A xx (i) = r)i6i t i, 
A yv (i) = Vifa.i and A xy (i) = j^flj. 

One may have already noticed that, the gauge trans- 
formation l|14|l is exactly the gauge transformation for 
the graviton if we view rank two tensor A^ as the lin- 
earized metric tensor on 2 dimensional space. Then the 
only gauge invariant ring exchange is the curvature ten- 
sor, which is R atl f3v Written in terms of linearized metric 
tensor, the curvature tensor is2£ 

Ranpv — l/2(A Q v,/n/3 + A^p^ av — A^ v ^ a p — A a fj^ v ). 

(15) 

Here the notation in general relativity has been used, 
A a v,jif3 — dfj,dpA au . 

In 2 dimensional space, although the curvature tensor 
has four indices, there is only one nonzero independent 
component, which is 

Rxyxy — ~^{^-xx.yy ~t~ Ayy XX 2^4.^^ ^^). (1^) 

At eighth order perturbation of t, a ring exchange term 
COS( K R X y X y ) is generated. Now the low energy effective 
Hamiltonian reads 

Heff = -tcos(R xyxy ) + —{E 2 XX + E^ + aE 2 X y) (17) 
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FIG. 5: The ring exchange of the 2 dimensional lattice. El- 
ementary hoppings are between site and nearest neighbor 
face center. Gauge invariant ring exchange term can be rep- 
resented as CQi,(Rxyxy\ it involves eight single hoppings of 
bosons. 



The cosine term in this Hamiltonian is ring exchange, 
t ~ t s /V 7 . If the original boson language is taken, one 
of the ring exchanges is shown in Fig|31 one can clearly 
see that the ring exchange term involves 8 independent 
nearest neighbor hoppings, this term only takes place at 
the eighth order perturbation of t. 

The first conclusion drawn from the gauge symme- 
try Ijl4|l is that, this system cannot be in the super- 
fluid phase. Because nonzero expectation value of the 
boson creation operator implies nonzero expectation of 
Aij . Any linear combination of Aij will break the local 
gauge symmetry, which is not allowed without matter 
fields. Another way to view this point is through the 
Hamiltonian. The hopping term Hi in i|10|) has global 
U(l) symmetry, which corresponds to the global conser- 
vation of total boson number. One might wonder whether 
this global 17(1) symmetry can be spontaneously bro- 
ken, i.e. the system becomes superfiuid. However, Hq 
opens a very big energy gap to excitations that do not 
conserve boson number. The big charge gap rules out 
the possibility of superfluidity. Because of the graviton 
gauge symmetry, the polarization of the collective excita- 
tion should automatically be the same as the gravitons. 
However, one crucial difference from the photon liquid 
phase is that, here the curvature tensor is the second or- 
der derivative of Aij , if there is a phase in which we can 
expand the cosine functions in (|17|) . i.e. there is a Gaus- 
sian phase or Gaussian fixed point, the dispersion relation 
of the gapless collective modes should be u> ~ k 2 . Un- 
fortunately, just like the monopole proliferation in 2+1 
dimensional QED, the topological defects in the current 
case also generally proliferate, and will gap out the gravi- 
ton excitations. 

The effect of the topological defects can be described 
in the dual formalism. Dual variable n and h can be 
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FIG. 6: The distribution of the dual variable h. h cannot all 
be integers over the whole 2 dimensional plane, instead, one 
quarter of h have to be half-integers. 



FIG. 7: The distribution of boson numbers on fee lattice. The 
link (in bold font) is shared by four plaquettes and two sites. 
Every site is occupied by three orbital levels, and on every 
plaquette there is one orbital level. 



defined as follows 

E xx = d y h, Eyy = d x h, 

E xy = —2d x d y h, 2R xyxy = n. (18) 

h and 7r are quantities defined on lattice i. One can check 
the commutator and see h and n are a pair of conjugate 
variables, i.e. [hi,iTj\ = iSij. Because the definition of 
the dual variable h only involves the second derivatives, 
the dual Lagrangian should be invariant under following 
transformation 

h->h + Ax + By + C (19) 

A, B, C are arbitrary constant integers. 
The dual Lagrangian now reads 

L du ai = (d T h) 2 - Pi h(d$ + d 4 + Aadld 2 y )h + ■■■ 

(20) 

The ellipses include vertex operators in this dual theory. 
Notice that, although E^ are all integers, the representa- 
tion of E xy in l|18l) contains factor 2 on the definition of 
the dual variables, thus h on some lattice sites have to be 
half integers, and the Berry phase will cause oscillation 
of the signs of the vertex operators on the lattice space. 
The distribution of h is shown in Fig. [B] The leading 
unoscillating vertex operators are 

Lvertea: = — « COs(47r/l) — ^ (cOs(Alt d x h) + C0s(47T<9j,/l)) . 

(21) 

This Lagrangian (|20|l looks exactly like the Lagrangian 
at the Rokhsar-Kivelson point for 2 dimensional dimer 
model on square lattic o 21 i 22 except for the vertex oper- 
ators. The leading kinetic term of this Lagrangian is 
proportional to fc 4 . The k 2 term is ruled out by the sym- 
metry (I19fl . Since k 2 term does not exist in the dual 
Lagrangian Q20|l . as long as 7 is tuned to zero and P4 is 
smaller than a critical value, the vertex operator is ir- 
relevant, and the system is a liquid phase without any 



order—. However, the vertex operator proportional to 7 
is relevant generically (and hence a term is dangerously 
irrelevant). The relevant vertex operator will drive the 
system into crystalline phase. The crystalline pattern 
close to the critical point with 7 = can be predicted 
from the field theory. The representations of density op- 
erators in terms of the low energy field variable are 

m,n 2 - 1 ~ (-If sin(27rcVi) + (-If sin(2ird y h), 
n-1 ~ -(-l) x sin(2TTd x h) - (-If sin(27n9 y /i), 
cos(R xyxy ) — (-l) y cos(2ird x h) + {-l) x cos{2-Kd y h). 

(22) 

When the vertex operators in l|21(l are relevant, the order 
parameters above will take nonzero expectation values. 
Thus the 7 operator tends to drive the system into the 
(it, 0) + (0, 7r) state for either particle density or plaquette 
density, depending on the sign of 7. Each case has 4 
degenerate ground states. 



IV. GRAVITON MODEL IN 3 DIMENSIONAL 
SPACE 

The model in 3 dimensional space is defined on the fee 
lattice, since physical quantities are defined on both sites 
and the centers of each plaquette. Let us assume there 
are 3 orbital levels on each site, and 1 orbital level on 
each face center. The particle number on the face center 
is denoted as n, and the particle numbers on sites are 
denoted as {ni 7 n 2 ,n 3 ) (FigO. 

The Hamiltonian for this system still contains three 
parts, H = Ho + Hi + H 2 - Hi is the nearest neighbor- 
hopping between sites and their nearest face centers, and 
also between adjacent face centers (notice that the adja- 
cent face centers have the same distance as the site and 
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its nearest face center). 

3 

Hi= £ £ -tblfij - g ififo + ft.c. (23) 

<j,j>a=l <i,i> 

, i/ 2 is the on site interaction H 2 = U(n—l) 2 , which fixes 
the average filling of the fee lattice. Hq is the interaction 
term involving links in all 3 directions. For example, for 
the link in (i,x), the interaction term reads 

Hq = V{n i+ i/2x + l/2z + n i+l/2x-l/2z + n i+l/2x+l/2y 
+ n i+l/2x-l/2y + 2711,,; + 2ni :i+x — 8) 2 . 

(24) 

. The links in y and z directions are treated similarly. If 
the bracket in equation 1(24(1 is expanded, it becomes the 
usual two body repulsion term. 

When Hq becomes the dominant term in the Hamil- 
tonian, it effectively imposes a constraint on the system. 
Again the best way to view this constraint is by introduc- 
ing a staggered sign and defining new variables, similar 
to the electric field in the dimer model discussed before. 
Let us define a rank-2 tensor Eij . The off-diagonal terms 
are defined on face centers as Eij = ?/ r (n— 1). n is located 
at one of the ij face centers; the diagonal term is defined 
on sites as En = r\ r [jii — 1) with £ = 1,2, 3. 77,- — ±1 and 
the distribution of sign rj r is shown in Fig. 0] 

After introducing the sign 77, the constraint effectively 
imposed by 1)24(1 can be compactly written as 

2d x E xx + d y E X y + d z E xz = 0, 

d X E X y + 2dyEyy + 8 Z Ey Z = 0, 

d x E xz + d y E yz + 2d z E zz = 0. (25) 

A violation of this constraint can be interpreted as a 
charged defect excitation and can drive the system into 
an ordered boson superfluid state after condensation. We 
will discuss the transition in section V. In the current 
section we only focus on the case when the constraint is 
strictly imposed. 

The constraint (|25|) requires the low energy Hamil- 
tonian to be invariant under the gauge transformation 
Aij — > Aij + difj + djfi. This is precisely the gauge sym- 
metry of the graviton in 3 dimensional space. Thus, the 
low energy physics can only involve the linearized curva- 
ture tensor. In the 3 dimensional space, the curvature 
tensor has 6 nonzero components, according to the sym- 
metry of the curvature tensor. Thus, now the effective 
low energy Hamiltonian reads 

Hring = £ —tl COS(Rijij) — £ t 2 COs(Rijik) 

ij,i^j ijk,i^j,j^=k,i^k 




(26) 

The cosine terms involving the curvature tensor are ring 
exchange terms generated by the nearest neighbor hop- 
ping. Ring exchanges in 1(26(1 are generated at the eighth 



order perturbation of the nearest neighbor hopping, i\, 
t 2 ~ t 8 /V 7 . All the lower order perturbations only gener- 
ate terms which do not comply with the constraint 125(1 . 
The ring exchange term cos(R xy xy) is the same as its 2 
dimensional counterpart, as shown in Fig. [5] 

In order to derive the correct Lagrangian, one needs 
to introduce a Lagrange multiplier Aio for the constraint 
((25(1 . The full Lagrangian reads 

L = y^(d T Ajj - diAjo - djA i0 ) 2 + ^ h cos(R i:jlj ) 

+ ^ hcos(Ri jik ) 

(27) 

The gauge symmetry can now be enlarged to quanti- 
ties defined in the 3+1 dimensional space-time: A^ v — > 
A^v + d^ v + d v f^ and ^4 o = 0, / = 0. In this sys- 
tem, the boson superfluid order is again ruled out by 
the gauge symmetry. Without crystalline order (proven 
later), the system is in a liquid phase with excitations 
which have the same gauge symmetry as the graviton. In 
the linearized Einstein gravity, after taking the traceless- 
transverse gauge, the spin-2 graviton has only two po- 
larizations. In our theory tracelessness was not imposed 
to A^ . Therefore there are three gapless collective exci- 
tations in the graviton phase, one of which is the scalar 
trace mode, the other two are described by traceless ma- 
trices, if the transverse gauge is taken, the two traceless 
modes exactly correspond to the two polarizations of the 
gravitons. 

Unlike the quantum dimer model, the curvature ten- 
sor is the second spatial derivative of A^ . If there is a 
Gaussian phase in which we can expand the cosines in 
equation ((27(1 . the gapless graviton mode in this Gaus- 
sian phase has a soft dispersion u> ~ k 2 . 

Whether the graviton excitations survive (or equiva- 
lently whether crystal order develops) can be studied in 
the dual theory. If we define the symmetric tensor £,j as 
£u = 2Eu, E.- L j — Eij,i 7^ j, the constraint i(25|) can be 
solved by defining the dual tensor hij as 

Eij — Ciabtjcddadchbd- (28) 

This is a double curl of the symmetric tensor hij. hij 
also lives on the sites and faces of this fee lattice. 

If checking carefully, one can notice that, the curvature 
tensor can also be written in the double curl form 

2Rxyxy — ^zab^zcd^a^cAbd-j 2R XZ xz — ^yab^ycd^a^cA^ci^ 
2Ry Z y Z — ^xab^xcd^a^cAbd-j 2R X yxz — ^yab^zcd^a^cA^ci^ 
2Ryxyz — €xab€zcd&a3cAbd 7 2R Z x Z y — ^xab^ycd^a^cA^d. 

(29) 

Therefore, this model is self-dual, as long as we define 
the dual variables hij in terms of £y = £iab£jcdd a d c hbd 
and its conjugate Try as follows: 

Rxyxy — ^zzi Ryzyz ^xxi 
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lx >xzxz yy i ^ 1Xj xzyz Jl xy> 

^Rxyxz — ^yzi ^-R-xyzy — ^xz- (30) 

According to the definition, -k^ is subject to the same 
constraint as Eij (|25|l . 

After introducing the dual variables h%j and nij, the 
dual Lagrangian reads 

Lduai = y^X®thjj — fyhjo - djhio) 2 - ^ Pi&ijij 
- P2R%ih H — 

(31) 

Rijki is the curvature tensor of hi j . hiQ is a Lagrange mul- 
tiplier, which is introduced for the constraint on iTij. The 
ellipses include possible vertex operators. Without the 
vertex operators, this theory is at a Gaussian fixed point 
and hence in an algebraic liquid phase with soft graviton 
excitations. If the vertex operators are relevant, they 
will destabilize the liquid phase and gap out the graviton 
excitation, and form crystalline order according to the 
Berry phase. The dual Lagrangian l|31|) is also invariant 
under the gauge transformation h^ v — > h^+d^f^+d^f^, 
with h 0Q — and f = 0. Therefore, any kind of vertex 
operator (for example cos(2Nwhij)) is not a gauge in- 
variant operator. The correlation function between two 
vertex operators at the Gaussian fixed point is zero or 
correlated at very short range (the correlation length is 
supposed to be roughly the inverse of the charge gap V in 
Hq (EU), so the Gaussian fixed point (also the algebraic 
spin liquid) is stable against weak perturbations of the 
vertex operators. Thus, this graviton phase is stable for 
the same reason as the photon phase, as discussed in the 
second section of this paper. 

V. PROPERTIES OF THE ALGEBRAIC LIQUID 
PHASE 

The most important feature of the algebraic spin liquid 
phase is the power law correlation between spin opera- 
tors. In the boson language used in this paper, it is the 
boson density operators which correlate algebraically. As 
discussed before, in the algebraic liquid phase, the boson 
density fluctuation around the average filling n can be 
written as (for instance) n\,i — n ~ r\it\ a b^\cdd a dchbd- 
The calculation of the equal time boson density corre- 
lation can be derived from the correlation functions be- 
tween hij. 

((m,i - n)(nij - »)) - ri i r) j d?d](h(i)h(j)) 

~r)iVjd 4 (-)~rii'nj^5- (32) 

r = \i — j\. Therefore the correlation functions between 
boson density operators fall off algebraically, with an ex- 



ponent higher than the exponent of the photon liquid 
phase. 

Because of the gauge symmetry, the operators have to 
be gauge invariant to have nonzero correlation functions. 
Therefore the correlation function between original boson 
creation operators b t are zero (or shortly correlated with 
a correlation length which roughly equals to the inverse of 
the charge gap V). The curvature tensors have nonzero 
correlations, and the correlators also fall off algebraically 
with the same exponent as that of the density correlator. 

The dynamics of the gapless liquid phase can be de- 
scribed by a new set of Maxwell's equations. Define the 
rank-2 tensor By = ^iabCjcddadcAbd, the dynamical equa- 
tions that describe this liquid phase can be derived di- 
rectly from the Lagrangian i|27|) , 

diSij 0, 

diBij = 0, 

8 t £ij - KtiabtjcddadcBbd = 0, 

d t Bij + Ke ia b£ 3 cddadc£bd = 0. (33) 

Charged excitations and topological defects are absent 
in these equations, thus the equations correspond to the 
Maxwell's equations in vacuum. If the constraint ll'">l) is 
softened, charge density and charge current have to be 
incorporated in equation l|33|) . 

The violation of constraint 125(1 corresponds to the de- 
fects, which carry gauge charges of gauge fields Aij. Be- 
cause the gauge field is a rank-2 tensor, the gauge charge 
should be vector field. The static vector charge field cou- 
ples to the gauge field in a similar way with the Gauss's 
law: diSij = —pj. For instance, if on one site, we increase 
n\ by one, it is equivalent to excite a pair of opposite 
gauge charges on the two x links sharing this site (Fig. 
|SJ); if on one plaquette center, is increased by one, 
gauge charges are created for four links sharing this pla- 
quette (Fig. |3 . Notice that since the constraint imposed 
by Hq in the original Hamiltonian (|24|l is an interaction 
between particles around each link, now the gauge charge 
field is defined on links. 

The defects should couple to the gauge field in a gauge 
invariant manner, in the rotor language, the coupling can 
be written as 

H e = -Y,t ^(8,6 f + d o ef ] - Aij ) + ■ ■ ■ (34) 

(e) 

6\ are the phase angles of gauge charge field creation 

operators. This coupling is gauge invariant since if 9^ 
is added by any arbitrary function /j, one can always 
eliminate this extra phase angle by gauge transformation 
— > + difj + djfi. In the Gaussian phase, the 

defects are gapped, and the phase angles 9^ are in the 
disordered phase. If the defects condense, i.e. the phase 
ang les 9<; e) take nonzero expectation values, the A^ will 
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FIG. 8: one of the defects carrying gauge charges. If m on 
one site is increased by one, it is equivalent to creating a pair 
of opposite gauge charges on two x links sharing this site. 
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FIG. 9: one of the defects carrying gauge charges, if on one 
plaquette center, rif is increased by one, gauge charges are 
created for four links sharing this plaquette. 

be gapped out through Higgs mechanism. However, the 
matter field phason modes will not be completely gapped 
out. This is because of the extra global total particle 
number conservation besides the gauge charge conser- 
vation. After the condensation of all the matter fields, 
the gauge charge conservation will be broken, and the 
gauge field is gapped out through the Higgs mechanism. 
Meanwhile, the spontaneous breaking of the global parti- 
cle number conservation guarantees the existence of the 
gapless excitation in the condensate, which is exactly the 
Goldstone mode. Notice that, although there are in total 
four flavors of matter fields (n%, 713 and n), different 
flavors of matter fields are mixed in the nearest hopping 
term Hi (I23[) . therefore only the total boson number is 
conserved. There should be only one Goldstone mode 
in the condensate. Similar situation takes place in the 
doped quantum dimer modeUS, where the holes carry 
both the gauge charge and the global t/(l) charge, hence 
the condensate of holes contains one gapless Goldstone 
mode. 

In the dual formalism, the topological defects can be 



introduced in the vertex operators of the dual formalism. 
The topological defects can be introduced in the dual 
vertex operators 

H m = - ^tcos{2TrN{d l 9 { m) + dj6l m) ) - 2nNh l:j ) 

(35) 

9^ are the phase angles of the creation operators of 
the topological defects. One can see the topological de- 
fects couple to the dual gauge potential hij in the same 
manner as the gauge charges couple to the original gauge 
potential Aij, i.e. the charge and topological charge are 
dual to each other, which is the same case as the duality 
between the electric charge and the magnetic monopolc. 
After the condensation of the topological defects 6^ , 
the gapless graviton excitation hij is gapped out by Higgs 
mechanism. However, in this situation, there is no extra 
global conservation of the topological defects, therefore 
in the condensate there will not be any gapless Gold- 
stone mode. The system is in the gapped Mott Insulator 
phase, with crystalline order determined by the Berry 
phase. Since our emphasis of this paper is about the sta- 
ble liquid phase, we will not get into the detailed analysis 
of the crystalline phase. 

After introducing the matter field and the topolog- 
ical defect, the semiclassical Maxwell's equations with 
charges are 

g.r.. - _J e ) 
diBij = -pf l \ 
d t £ i3 - Ke iab e jcd d a d c Bbd = + J^ 3 , 
d t Bi 3 + KeiabtjcddadcSbd = Jj + J- m ' j . (36) 

(e)i 

The current ■/_• represents the current with the zth 
component of matter field flowing in j direction. The 
conservation law of the gauge charges, as well as the topo- 
logical defects reads 

dt pf = -d, Jp - d t Jj e) \ 
%>J m) = -9i4 mW - di Jj ro)i . (37) 

One can clearly see that, because of the flavor mixing 
in Hi, the vector charge density is not individually con- 
served. The derivative diJ^ 1 plays the role of torque 
density of the vector charge density. This is similar to 
the spin currents^, which is also a tensor current, and the 
conservation equation of the spin current involves torque 
densitySi, due to the fact that spin is not conserved in a 
system with spin-orbit coupling. 
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VI. CONCLUSIONS, EXPERIMENTAL 
REALIZATIONS AND EXTENSIONS 

In this work we constructed models which give rise to 
collective excitations analogous to the gravitons. In 2+1 
dimensional system, the graviton excitations are unstable 
against vertex operators in the dual formalism; in 3+1 
dimensional space, there is a stable algebraic boson liquid 
phase which contains gapless graviton excitations with 
soft dispersion. The graviton liquid phase is self-dual, 
and the stability of the algebraic phase is guaranteed by 
the large gauge symmetry on both sides of the duality. 
The dynamics of the phase is described by a new set of 
Maxwell's equations. 

We proved the algebraic phase is stable, i.e. the re- 
alization of this phase does not require any fine-tuning. 
However, the original bosonic model is of a special form, 
its precise realization in experimental systems requires 
more efforts. Here we only consider the possibility of its 
realization in the cold atom system trapped in an optical 
lattice. The original boson model contains standard near- 
est neighbor hopping between particles, and also contains 
two body interactions, both on-site and off-site. The on- 
site repulsion can be obtained from the s-wave scattering 
between bosons on the same site, and the off-site repul- 
sion has been shown recently to exist in Chromium atom 
condensate due to the long range dipole interactions^*^. 

The interaction between the three orbital levels on 
each site and the bosons in the plaquette centers is very 
anisotropic (|24|) . this could be due to the anisotropy of 
the orbital level spatial wave functions. In the transi- 
tion metal oxides materials, the anisotropy of the t2 g 
level electron wave functions gives rise to the Khali- 
ullin model2L2& which was supposed to explain the or- 
bital liquid phenomenon, in which the orbital moment is 
quenched by quantum fluctuation — . Recently it has been 
proposed that, although the ground state of each site of 
the optical lattice is s-wave, if the particles are pumped to 
the first excited three fold degenerate p-wave states, the 
particles will maintain in the excited states for a consid- 
erable time, long enough for the particles to equilibrate^, 
therefore the p-wave particles will first form equilibrium 
state before they drop to the s-wave ground states. Each 
of the three fold degenerate p-level wave functions only 
extends in one direction of the 3 dimensional space. Since 
the wave function overlap is anisotropic in space, the par- 
ticle on each orbital level interacts more strongly in one 
certain direction, as long as the optical trap on each site 
is not too deep, i.e. the wave function is not too localized 
on each site. The anisotropic interaction is required in 
our problem. The remained problem of the experimen- 
tal realization is how to realize the fee lattice with laser 
beams, and how to pump particles on all the sites of the 
fee lattice to the excited p-wave states. 

So far the models we have considered are purely 
bosonic models. After the search for unconventional 
bosonic phases, one might wonder if non-fermi liquid can 
be obtained in a similar way for fermionic systems. There 



has been a great deal of study on treating the constraint 
of no double occupancy of Hubbard model by introducing 
[7(1) gauge field, and the bosonic holons and fermionic 
spinons are interacting with this U(l) gauge field. The 
interacting system has been proposed to explain the non- 
fermi liquid behavior of the normal state of High Tc cu- 
perates (for instance, referenced). In this kind of theo- 
ries, the emergent U (1) gauge field plays the most crucial 
role. For instance, due to the scattering from the gauge 
fields, the scaling of the resistivity significantly deviates 
from the T 2 law of normal fermi liquid. The emergence 
of this gauge field is exactly due to the local constraint 

In our current work, if all the particles in the origi- 
nal model are fermions, the constraint imposed by Hq 
can also be solved by introducing rank-2 tensor gauge 
field Aij, and when the gauge field is in its deconfined 
phase, the system can be viewed as fermions interacting 
with gapless soft graviton mode . The behavior of the 
fermions is expected to be non-fermi liquid. 

In the whole paper, our goal has been carefully lim- 
ited to the discovery of a new type of algebraic spin liq- 
uid. However, one can consider the issue of "quantum 
gravity" as an extension of this work. Indeed, if in a 
theory the gauge symmetry of the graviton can emerge 
at low energy physics, the theory might be a candidate 
of a new possibility of quantum gravity. Several other 
systems in condensed matter physics have been proposed 
to be related to gravity, for instance, spin-2 particles are 
supposed to exist at the edge states of 4 dimensional 
quantum Hall models. However, in these systems the 
gauge symmetry (which is very crucial for gravitons) was 
not an emergent property. In our work, the collective ex- 
citations automatically have the gauge symmetry of the 
graviton. However, the dispersion relation is quadratic, 
this is due to the fact that the gauge invariant operator, 
the curvature tensor is the second derivative of Aij . The 
ring exchange term generated from the nearest neighbor 
hoppings is cos(i?), after the expansion, the leading term 
is R 2 , which is proportional to fc 4 . 

We can make the graviton dispersion linear by intro- 
ducing by hand a quasi-gauge invariant term to the low 
energy Lagrangian : 

cddad c A bd . (38) 

This term is not completely gauge invariant, instead, it 
is gauge invariant up to a boundary term. This is very 
analogous to the Chern-Simons field theoy of the Quan- 
tum Hall effect^. However, this term cannot be gener- 
ated from the microscopic boson model with only nearest 
neighbor hopping. It is expected that, by coupling to the 
matter fields, and the matter fields form certain special 
state, the k 2 term (|38|l can be generated from integrat- 
ing over the matter fields in the partition function, just 
like how the Chern-Simons field theory is obtained in the 
Quantum Hall state of electrons. Recently, some other 
authors have proposed a bosonic model which is simi- 
lar to ours, and the k 2 term (|38H is claimed to exist in 
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the long scale physics, therefore a graviton with linear dispersion is supposed to exists*. 



1 P. W. Anderson, Science 235, 1196 (1987). 

2 G. Baskaran, Z. Zou, and P. W. Anderson, Sol. St. Comm 
63, 973 (1987). 

3 E. Fradkin and L. Susskind, Phys. Rev. D. 17, 2637 (1978). 

4 R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 
(2000). 

5 L. Balents, M. P. A. Fisher, and S. M. Girvin, Phys. Rev. 
B. 65, 224412 (2002). 

6 D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 
2376 (1988). 

7 A. Vishwanath, L. Balents, and T. Senthil, Phys. Rev. B 
69, 224416 (2004). 

8 A. M. Polyakov, Nucl. Phys. B. 120, 429 (1977). 

9 A. M. Polyakov, Gauge Fields and Strings, Harwood Aca- 
demic, New York (1987). 

10 R. Moessner and S. L. Sondhi, Phys. Rev. B. 68, 184512 
(2003). 

11 X. G. Wen, Phys. Rev. B. 68, 115413 (2003). 

12 M. Hermele, M. P. A. Fisher, and L. Balents, Phys. Rev. 
B. 69, 064404 (2003). 

13 T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and 
M. P. A. Fisher, Science 303, 1409 (2004). 

14 T. Senthil, L. Balents, S. Sachdev, A. Vishwanath, and 
M. P. A. Fisher, Phys. Rev. B. 70, 144407 (2004). 

15 M. Hermele, T. Senthil, M. P. A. Fisher, and a. et, Phys. 
Rev. B. 70, 214437 (2004). 

16 J. Goldstone, A. Salam, and S. Weinberg, Phys. Rev. 127, 
965 (1962). 

17 C. Xu, Condmat/0602443 (2006). 

18 L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and 
K. Sengupta, Phys. Rev. B 71, 144508 (2005). 



1 L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and 
K. Sengupta, Phys. Rev. B 71, 144509 (2005). 

20 C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravita- 
tion, WH Freeman Press, San Francisco (1973). 

21 C. L. Henley, J. Stat 89, 483 (1997). 

22 E. Fradkin, D. A. Huse, R. Moessner, V. Oganesyan, and 
S.L. Sondhi, Phys. Rev. B. 69, 224415 (2004). 

23 S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301, 
1348 (2003). 

24 D. Culcer, J. Sinova, N. A. Sinitsyn, T. Jungwirth, A. H. 
MacDonald, and Q. Niu, Phys. Rev. Lett 93, 046602 

(2004) . 

25 A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and 
T. Pfau, Phys. Rev. Lett 94, 160401 (2005). 

26 J. Stuhler, A. Griesmaier, T. Koch, M. Fattori, T. Pfau, 
S. Giovanazzi, P. Pedri, and L. Santos, Phys. Rev. Lett 95, 
150406 (2005). 

27 G. Khaliullin and S. Maekawa, Phys. Rev. Lett 85, 3950 
(2000). 

28 G. Khaliullin and S. Okamoto, Phys. Rev. Lett 89, 167201 
(2002). 

29 B. Keimer and et al., Phys. Rev. Lett 85, 3946 (2000). 

30 A.Isacsson and S. M. Girvin, Phys. Rev. A 72, 053604 

(2005) . 

31 P. A. Lee and N. Nagaosa, Phys. Rev. B 46, 5621 (1993). 

32 S.-C. Zhang and J. Hu, Science 294, 823 (2001). 

33 S. C. Zhang, T. H. Hansson, and S. Kivelson, Phys. Rev. 
Lett 66, 82 (1989). 

34 Z.-C. Gu and X.-G. Wen, gr-qc/0606100 (2006). 



